Preparation of Subradiant States using Local Qubit Control in Circuit QED 
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Transitions between quantum states by photon absorption or emission are intimately related to 
symmetries of the system which lead to selection rules and the formation of dark states. In a circuit 
quantum electrodynamics setup, in which two resonant superconducting qubits are coupled through 
an on-chip cavity and driven via the common cavity held, one single-excitation state remains dark. 
Here, we demonstrate that this dark state can be excited using local phase control of individual 
qubit drives to change the symmetry of the driving held. We observe that the dark state decay 
via spontaneous emission into the cavity is suppressed, a characteristic signature of subradiance. 
This local control technique could be used to prepare and study highly correlated quantum states 
of cavity-coupled qubits. 
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Symmetry properties of a quantum system interacting 
with a radiation field provide information about possi- 
ble transitions within the system. Symmetry operations 
such as translation, rotation or reflection which leave 
the system invariant, lead to selection rules in molecular 
and solid-state systems pQ. For an ensemble of identi- 
cal atoms, the symmetry under permutation of particles 
allows only for transitions between symmetric collective 
states [2j [3]. Such highly-entangled Dicke-states, like the 
single-excitation W-state, have attracted a lot of atten- 
tion in the field of quantum information processing with 
trapped ions [4 , optical photons [5] or superconducting 
qubits [6 due to their robustness under decoherence [7] 
and particle loss [8j[9]. Moreover, collective states have 
been used for quantum information storage in atomic 
memories [10 J. The particular symmetry of Dicke states 
affects also their decay. Spontaneous emission can be en- 
hanced for superradiant states [2J[3], or inhibited, if the 
symmetry of the state does not allow for the emission 
of a photon. This effect is known as subradiance and 
is closely related to the concept of decoherence-free sub- 
spaces, regions in Hilbert space which are not affected by 
decoherence and therefore appealing for quantum infor- 
mation processing [11]. Though theoretically well stud- 
ied p21 D3] , subradiant states are difficult to realize and 
experimental evidence of subradiance is rare Q~4j [15] . 

Here we present a method to prepare two qubits in an- 
tisymmetric Dicke states and demonstrate their subradi- 
ance in circuit quantum electrodynamics (QED) [T6| ITT] . 
In this architecture, superconducting artificial atoms are 
coupled to a common field mode of a planar microwave 
cavity. Strong resonant coupling of individual qubits [TSl - 
[20] and qubit-ensembles [21] to single microwave photons 
has been achieved. Moreover, cavity-mediated interac- 
tions between distant qubits [22-24] form the basis for 
on-chip quantum information processing [25] [26] using 
entangled states of currently up to three qubits [27] . In 
many experiments, the qubits are detuned from the cav- 
ity resonance. In this dispersive regime, the cavity can 
be used for the preparation and read-out of entangled 



states [28, 29 , and cavity induced radiative decay of the 
qubits due to the Purcell effect [30 is reduced. This 
enhanced spontaneous emission can be modified by the 
qubit-cavity detuning [31] and is the dominant relaxation 
mechanism close to resonance. We have observed subra- 
diance by preparing an antisymmetric two-qubit state in 
this regime using local phase control to change the sym- 
metry of the driving field [Fig. [TJa)]. 

We consider two qubits resonant with each other but 
not with the cavity, modeled by a generalized Tavis- 
Cummings Hamiltonian [32] . 

H^c/h = oj r a)a + uo q J z + g(aJ+ + a) J_). (1) 

where the dominant cavity mode is at frequency uo r and 
the qubits are at frequency uj q . The operators J z = 



n U) 



1/2 ^ /i a z and J± = J^ <j± are collective spin op- 



erators [3] with a± = (ax ± ia y z) )/2 and Pauli opera- 
tors (Tx,y,z for the individual qubits (i = 1,2). a^) is the 
annihilation (creation) operator of the field interacting 
with the qubits with equal coupling strength g. For the 
single-excitation manifold the eigenstates of H^c are 

|Va> = |0;V->, ( 2 ) 

|Vv) = cos6 m \l;gg) + sin# m |0; ^+), 
\ip s ) =sin0 m |l;£flr) - cos# m |0; ^+), 

where the mixing angle 6 m is given by cos2# m = 



— A/a/4(v / 2^) 2 + A 2 . |n; ) is a state with n photons 

in the resonator, and \ip±) = (\ge) ± \eg))/y/2 are the 
symmetric and antisymmetric Bell- states. The qubit- 
resonator detuning A = uo q — uo r is chosen negative. 
The antisymmetric state \ip a ) - comprising the antisym- 
metric qubit state |^_) - does not couple to the cavity 
field. Only the symmetric qubit state |^ + ) mixes with the 
field to form hybridized matter- field states \ip r ) and \i(j s ) 
[Fig. [TJb)] with enhanced collective coupling strength 
\ /f Ng forN = 2 [21] . Generally, for an 7V-qubit ensemble 
there are always two hybridized bright states and N — 1 
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FIG. 1: (a) Schematic drawing of a cavity QED setup with in- 
dividual phase control of the driving field for each (artificial) 
atom, (b) Energy level diagram of the two-qubit system cou- 
pled dispersively to a common cavity field. Symmetric states 
are indicated by thick solid (green) lines and the antisymmet- 
ric state is identified by the thick dashed (red) line. 



uncoupled, dark states. Only qubit states which reflect 
the symmetry of the cavity mode (i. e. belong to compat- 
ible irreducible representations of the permutation group 
[T2] ) couple to the single-photon field. In our case, the 
coupling to the first harmonic cavity mode has the same 
sign for both qubits (g^ ~ g^) [24]. Therefore, the 
symmetric qubit state 1^+)? which is invariant under per- 
mutation of the qubits, couples to the field. 

For our experiments, two superconducting transmon 
qubits [33] have been integrated into a coplanar nio- 
bium resonator on a sapphire substrate [Fig. |2ja,b)]. 
The qubits have similar Josephson energies Ej/h ~ 
37.6 GHz, charging energies Ec/h « 285 MHz and cou- 
pling strengths g/2ir ^116 MHz to the first harmonic 
mode of the microwave transmission line resonator. The 
resonator frequency is uj r /27r = 6.937 GHz and its decay 
rate is k/2tt = 3.01 MHz. In the dispersive regime where 
8 m « 7r, the photonic contribution |1; gg) to the symmet- 
ric state \i/j 3 ) is small (of order y/2g/A) and the state has 
pre-dominantly qubit character. It is however, shifted in 
energy by 2 J = 2g 2 /A corresponding to the dispersive 
J-coupling discussed in Refs. [22] [24] [Fig. [ljb)]. The 
antisymmetric wave- function \ip a ) has no photonic com- 
ponent and thus experiences no Lamb-shift of its energy. 
In this description the qubit-qubit coupling J can be un- 
derstood as the collective Lamb-shift (\/2g) 2 / 'A = 2J of 
the symmetric state \ip s ). 

The symmetry of these states is also reflected in se- 
lection rules for electric dipole transitions. In fact, for a 
drive applied directly to the cavity only, transitions from 
the ground to the symmetric bright state \i(j s ) are al- 
lowed, while transitions to the antisymmetric dark state 
\ip a ) are forbidden [24]. Again, this is due to identical 
microwave fields at the qubit positions for a drive close 
to the first harmonic cavity mode, which conserves the 
symmetry under permutation of qubits. Then only tran- 
sitions within the class of symmetric states [Fig. [lib); 
solid green lines] are allowed [12] . This constraint can be 
overcome by addressing the qubits individually via capac- 



itively coupled charge lines [28] and tuning the relative 
phase (j) of the microwave drive at the qubit positions. 
When choosing a relative phase of <p = 7r, the opposite 
sign of the local fields results in allowed transitions to the 
antisymmetric state [42] . 

In the rotating frame, the drive acting on the individual 
qubits with frequency ujd and coupling e is 



H d = he{a^+^a^) 



h.c. 



(3) 



where £ is the amplitude imbalance at the individual 
qubits. Starting in the ground state, the drive Ha can 
induce transitions to the state |^), if the matrix ele- 
ment Q(ip) = \(0;gg\H d \ilj}\/h is non-zero. ^(^ dark ) = 0, 
therefore, defines the dark state for given drive imbalance 
£ and relative phase (j). In the dispersive regime, the ma- 
trix element for the symmetric and antisymmetric state 



fi(Vy a ) = ev /(l + £ 2 ±2£cos0)/2. 



(4) 



For equal drive amplitudes (£ = 1) and zero relative 
phase ((f) = 0), which corresponds to a drive applied 
to the cavity in the vicinity of the first harmonic mode, 
£l(ip a ) vanishes and the antisymmetric state \ip a ) remains 
dark [24 J. For (j) = tt and £ = 1 however, Cl(ip a ) is max- 
imal, while the transition rate fl(ip 3 ) to the symmetric 
state vanishes. The transition can thus be enabled, or 
disabled, by adjusting the relative phase appropriately. 

We have spectroscopically measured the transition am- 
plitude as a function of relative phase <\> in the vicinity 
of the bare qubit frequency uj q /27T = 6.647 GHz. To 
control (/>, the local microwave fields are generated by 
a single microwave source operating at the carrier fre- 
quency cjlo = ^d+^iF- We use two in-phase/quadrature 
(IQ) mixers to generate sidebands of the carrier signal 
at the frequency uod- The signal (IF) at the interme- 
diate frequency cjif = 150 MHz is synthesized with an 
arbitrary waveform generator and applied to each mixer 
with relative phase <\> [Fig. [2ja)]. The phase (j) of the lo- 
cal qubit drive fields can then be controlled with high 



^"^^ ADC 




@i Q(g )--rTHv> 



FIG. 2: Setup and micrograph of the sample with a transmis- 
sion line resonator (green) and two transmon qubits (blue) 
addressable via local charge lines (yellow). 
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FIG. 3: Spectroscopy of dark and bright state as a function 
of the relative phase <j> between drives, (a) Normalized inte- 
grated transmission amplitudes S n = S/Sm ax of the antisym- 
metric state ipa (black solid dots) and the symmetric state ip s 
(red open diamonds) at the frequencies cj a /27r = 6.647 GHz 
and (jO s /2tv = 6.578 GHz, respectively, as indicated by the ar- 
rows in (b). S n is normalized to unity at the maximum value 
of both curves S'max- The lines indicate fits to the expected 
steady-state populations, (b) Transmission amplitude as a 
function of frequency va and relative phase 0. 



precision by the phase of the IF-signals. Pulsed spec- 
troscopy is used [34] [35] , where the measurement tone is 
switched on after a 500 ns saturation pulse [Fig. |3ja)]. 
The difference between the transmitted signal s(t) and 
the signal with both qubits in their ground state s gg (t) 
integrated over t m = 520 ns yields the transmission am- 
plitude S = j m {s(t) — s gg (t)}dt. S corresponds to a 
measurement of the steady-state population of the single- 
excited qubit states when appropriately normalized [29] . 
Two spectroscopic lines, varying in amplitude as a 
function of relative phase </>, are observed at cj a /27r = 
6.647 GHz and uj s /2tt = 6.578 GHz. These corre- 
spond to the states \ip a ) and \i/) 8 }, respectively. As 
expected, the zeros of the populations of the symmet- 
ric and the antisymmetric state are 180° out-of-phase 
[Fig. [5]. The measured transmission amplitudes are in 
good agreement with the steady-state population of the 
excited states calculated using the dissipative Bloch equa- 
tions for a two- level system driven at zero detuning [36] , 
S = S {1 - 1/(1 + TiT 2 ft(V0 2 )}/2. The scaling factor 
So takes account of the integrated voltage at the analog- 
digital converter. The dephasing times T2 jS = 322 ±26 ns 
and T 2 , a = 526 ± 36 ns of \ifj a ) and \ip s ) are determined 
in independent Ramsey- fringe experiments. The energy 
relaxation times T\ i3 and Ti >a are measured in separate 
time-delay experiments and are discussed below. From 
independent fits of S to both curves we obtain the rel- 
ative phase difference between the zeros of the popu- 



lations of the symmetric and the antisymmetric state 
</> s - ^a = 7T - 0.3 ± 0.03 rad using Q(ip 3 / a ) from Eq. (Eh. 
The deviation of the measured phase difference from it 
is simply caused by the difference in cable lengths of the 
two qubit drive lines. 

Using this method, it is possible to verify the subra- 
diant character of the antisymmetric state by testing its 
resilience to cavity-induced Purcell decay [30 , which is 
caused by the indirect coupling of the qubits to the envi- 
ronment via the cavity. According to Fermi's golden rule, 
the induced voltage fluctuations ex (a^ + a) of the cav- 
ity field lead to a decay rate 7^ = «|(0;^|a|^)| 2 to the 
ground state p~8j [33] . The total decay rate is then given 
by 7 = 7z + 7k with the intrinsic, non-radiative decay 
rate 7^. Although the Purcell decay can be made small 
by operating the qubits in the dispersive regime, where 
Ik ~ (9/ 'A) 2 ft, or by using advanced circuit designs [37] , 
it cannot be fully avoided for single transmon qubits. For 
the dark state however, the matrix element |(0; p^lal^)! 
vanishes completely, since by symmetry \ip a ) has no pho- 
ton admixture. In other words, destructive interference of 
the photons emitted from either qubit leads to a suppres- 
sion of the spontaneous emission process and the dark 
state is protected against Purcell decay. Note that a two- 
island transmon design, in essence the integrated version 
of the two-qubit/cavity circuit used in the experiments 
discussed here, also provides Pur cell-protect ion based on 
the formation of an intrinsic dark state [38l EH] . 

In order to observe subradiant Purcell-protection in 
a regime where radiative losses dominate over intrinsic 
qubit losses (7^ > 7^), we have detuned the qubits from 
the first harmonic mode by A/2tt = 290 MHz ~ 2.5g. 
At this frequency, the lifetime of the symmetric and an- 
tisymmetric state, as well as the decay rates of the in- 
dividual qubits have been measured. A delayed mea- 
surement pulse technique has been employed, where we 
apply a ir pulse resonant with the respective transition 
frequency and delay the time At before applying the 
read-out pulse. The lifetimes of single qubit excitations 



Ti 



ge 



401 ± 16 ns and Ti 



eg 



364 ± 16 ns at this 



frequency are comparable to the bright state lifetime 
of Ti 5 g = 368 ± 30 ns. In our measurements, the ef- 
fect of superradiant decay is masked by the intrinsic de- 
cay rate and by dephasing acting locally on individual 
qubits. In contrast, the measured dark state lifetime 
Ti, a = 751 ± 13 ns exceeds these values by a factor of 
two - a clear signature of subradiance that demonstrates 
the decoupling of the antisymmetric state from the cavity 
environment and, as a consequence, its enhanced stabil- 
ity. The population decay versus time of both the bright 
and the dark state is plotted in Fig. pla). 

The lifetime of the dark state is shown at different de- 
tunings, along with the lifetimes of the bright state and 
the uncoupled single qubit states in Fig. [4fb). It is veri- 
fied in numerical master equation simulations of the dis- 
sipative dynamics for ip a [Fig. Efb); solid black line] and 
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FIG. 4: (a) Population of the antisymmetric dark (ip a , black 
dots) and the symmetric bright (ip s , red diamonds) state ver- 
sus time at A/2tt = 290 MHz (gray line in (b)) with exponen- 
tial fits to the data. The inset shows the energy relaxation 
from ipa and ip s to the ground state, where only the symmet- 
ric state ip s is affected by the Purcell decay j^ . (b) Measured 
decay times of ip s (red diamonds) , ip a (black dots) and the un- 
coupled qubit states eg (orange points) and ge (green crosses) 
as a function of detuning A. Exponential fits to numerically 
simulated populations are shown for ip s (dashed line) and ip a 
(solid line). 



ip s [Fig. plb) ; dashed red line] that the decreasing lifetime 
of the dark state at small detunings is caused by local de- 
phasing. For these simulation, a pure dephasing time of 
^2,0 = 880 ns has been determined in Ramsey- fringe ex- 
periments of the uncoupled qubits with an intrinsic decay 
time of 1.37 /is obtained from a fit to Ti jeg (A). 

The enhanced dark-state lifetime can be used for quan- 
tum computation. In fact, the logical qubit formed by 
the ground and the dark state realizes a decoherence- 
free subspace, which is insensitive to cavity dissipation 
as well as to global dephasing acting simultaneously on 
both qubits. Note however, that the subspace spanned 
by |0; gg), \ip a ) and the doubly excited state |0; ee) forms 
a weakly nonlinear qubit with anharmonicity 2 J, which 
limits the shortest preparation time without pulse opti- 
mization jiH] to ~ 1/(2J). 

More generally, this local control technique may allow 
us to excite highly entangled Dicke states belonging to 
different symmetry classes with a single microwave pulse 
conditioned on the choice of phases between individual 
drives. Moreover, the possibility to address states of dif- 
ferent symmetry classes of multi-qubit systems can be 
used to encode information in collective qubit states. For 
readout, they can be transformed into entangled states 
in the computational basis by rapidly detuning the qubit 
transition frequencies. 

In conclusion, we have demonstrated a method to pop- 
ulate dark states of a two-qubit system in a circuit QED 
setup. The transitions to either dark or bright two-qubit 
states can be selected by adjusting the relative phase be- 
tween individual qubit drives, thus changing the sym- 
metry of the field and enforcing a symmetry-induced se- 



lection rule. We apply this technique to demonstrate 
Purcell-protection of the subradiant dark state against 
spontaneous emission. An extension to more qubits could 
provide further insight into the unitary and dissipative 
dynamics of multi-particle quantum states that can be 
directly prepared in the coupled qubit basis. Controlling 
the symmetry of the radiation field is, therefore, a vi- 
able method for preparation of states that are otherwise 
difficult to realize. 
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